This dataset accompanies Cheng, Hensley, Chang, and Doré (2025)
(arXiv: 2411.12801; published in ApJ).  

We present a new dust model for the North and South Galactic Caps (NGC/SGC),
constructed using gas tracers across multiple phases. 
Our approach employs a data-driven clustering algorithm to derive 
a new set of template bases separately for NGC and SGC. 
Specifically, our model consists of two H I component maps 
(LVC- and LVC in NGC; LVC- and LVC+ in SGC) and one H II map.  

This dataset includes maps of each template component, 
as well as model maps obtained by fitting our templates to the E(B-V) map from SFD
and to the three high-frequency Planck maps (353, 545, and 857 GHz). 
Residual maps from the Planck map fits are also provided 
(bottom panel of Figure 6 in the paper).  

The maps are provided in Galactic coordinates in HEALPix format 
with RING ordering and Nside = 1024.  

Below is an example script to access and visualize the 857 GHz Planck residual map,
corresponding to the bottom-right panel of Figure 6 in the paper.  

==============
from astropy.io import fits  
import healpy as hp  

hdul = fits.open('Cheng25_DR_maps.fits')  
hdul.info()  
data = hdul[1].data  
residual_map = data['Planck857_residuals']  
hp.orthview(residual_map, rot=[0,90,0], cmap='coolwarm')  
==============

Cheng25_DR_maps.fits contains the following columns:  

'NGC_mask': Mask for NGC  
'NGC_HI_IVC-' [cm^-2]: H I column density map for the NGC IVC- component  
'NGC_HI_LVC' [cm^-2]: H I column density map for the NGC LVC component  
'NGC_HII' [cm^-2]: H II map for NGC  
'SGC_mask': Mask for SGC  
'SGC_HI_LVC-' [cm^-2]: H I column density map for the SGC LVC- component  
'SGC_HI_LVC+' [cm^-2]: H I column density map for the SGC LVC+ component  
'SGC_HII' [cm^-2]: H II map for SGC  
'EBV_SFD' [mag]: All-sky SFD E(B-V) map at Nside = 1024, scaled by 0.884  
'EBV_model' [mag]: E(B-V) model map from fitting our templates to the SFD E(B-V) map  
'Planck353' [MJy/sr]: Planck 353 GHz map  
'Planck545' [MJy/sr]: Planck 545 GHz map  
'Planck857' [MJy/sr]: Planck 857 GHz map  
'Planck353_model' [MJy/sr]: Planck 353 GHz model map from fitting our templates to the Planck 353 GHz map  
'Planck545_model' [MJy/sr]: Planck 545 GHz model map from fitting our templates to the Planck 545 GHz map  
'Planck857_model' [MJy/sr]: Planck 857 GHz model map from fitting our templates to the Planck 857 GHz map  
'Planck353_residuals' [MJy/sr]: 353 GHz Planck map - 353 GHz Planck model map  
'Planck545_residuals' [MJy/sr]: 545 GHz Planck map - 545 GHz Planck model map  
'Planck857_residuals' [MJy/sr]: 857 GHz Planck map - 857 GHz Planck model map